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1 .  INTRODUCTION 


k 


•< 


Let  X  =  (x  ,...,x  )  be  a  random  vector  with  joint  d.f.  (distribution 

1  p  J 

function)  F,  i-th  marginal  d.f.  F.,  (i,j)-th  marginal  d.f.  Fj .  and  i-th 
marginal  density  function  f^.  We  denote  the  i-th  marginal  quantile  function 
by 

£^(q)  =  F^(q)  =  inf  {xtF^x^q}  ,0  <  q  <  1  (1.1) 

and,-  for  convenience,  a  specific  quantile  say  the  q^-th  of  F.  by 

=  5.(q.).  (1.2) 

Further . let 


rlij(q»f)  =  Fij(Ci(q).Cj(r)) 


(1.3) 


and  denote  for  given  q^  and  q^ 


°ij "  ni J  -  Vj '  hj'-VV'V'r 

The  parameters  (1.1)  -  (1.4)  defined  above  refer  to  the  d.f.  of  X. 
Now  let 


(1.4) 


X,  =  (x 


li’  •  ‘ 


•V’ 


i  =  1 


(1.3) 


be  n  independent  copies  of  X  and  denote  the  empirical  d.f.  of  {X^ , i=l , . . . ,n} 
by  F^n^  and  the  corresponding  i-th  and  (i,j)-th  marginal  distributions  by 
F^n^  and  F^  ™  ^  respectively.  We  denote  the  quantities  (1.1)  -  (1.4)  defined 
in  terms  of  F(n),  F^n)  and  F^  by 

,  (  n  ) .  .  _  (  n  )  .  (  n  ) 

(q),  9^  and  a 


(1.6) 


2 


^(q),  0  <  q  <  1,  i  =  1 . .  (1.9) 

The  asymptotic  distributions  of  the  empirical  quantle  process  (CsorgH  and 
Revesz  (1981))  and  of  a  fixed  set  of  specified  quantiles  (Mosteller  (1946)) 
in  one  dimension  are  well  known. 

Of  particular  interest  is  the  joint  asymptotic  distribution  of  the 
marginal  sample  medians 

(Cjd) . s  (D)  d.io) 

using  which  we  develop  tests  of  significance  for  the  population  medians 
analogous  to  tests  for  the  means  in  the  multivariate  case  (see  Rao  (1973), 
pp.  543-573).  An  early  work  on  the  joint  asymptotic  distribution  of  the 
sample  medians  is  due  to  Mood  (1941),  where  he  assumes  the  existence  of  the 
density  function  for  the  vector  variable  X.  We  obtain  the  distribution  in 
the  general  case  in  a  form  convenient  for  practical  applications. 


3 


2.  DISTRIBUTION  OF  THE  MARGINAL  SAMPLE  QUANTILES 

We  prove  the  following  theorem  concerning  the  joint  asymptotic  distribution 
of 

(0J.....0  )  =  UjCqj) . 1  (qp))  (2.1) 

the  sample  q^-th, . . . ,q  -th  quantiles  of  the  marginal  empirical  distributions  of 
x^,...,x  respectively. 

Theorem  2.1.  Let  be  continuously  twice  dif ferent iabie  iti  a  neighborhood 
of  0^  and  6^  =  f  ( T.-(q j )  )=f A 6 .  )  >  0,i=l,...,p,  where  f.  denotes  the  derivative  of 
F^.  Then  the  asymptotic  distribution  of 

y  '  =  . 0  -0  )  (2.2) 

n  11  P  P 

is  p-variate  normal  with  mean  vector  zero  and  variance-covariance  matrix 


/ 


qjU-Hj)  0 


12 


JjL 


Z  = 


6162 


6,6 

1  P 


(2.3) 


\ 


Pi. 


_£2_ 


6  6, 

P  1 


6  6» 
P  2 


q  ( 1 -q  ) 
_E _ P_ 


where  are  as  defined  in  (1.4). 


Proof .  By  Bahadur's  representation  of  the  sample  quantiles  (see  Bahadur 
(1966)). 


(log  n)  *n^^  |  (9j  -0j  )  -  <S11(r1~q^)|  — ^  0,  i  =  l 


(2.4) 


where  r^  =  F^n^(0^).  Then,  it  follows  that 


have  the  same  asymptotic  distribution.  By  the  multivariate  central  limit 


theorem,  z^  weakly  converges  to  a  p-variate  normal  distribution  with  mean  vector 


zero  and  covariance  matrix  as  given  in  (2.3).  This  proves  Theorem  2.1. 


For  practicle  applications  we  need  a  consistent  estimate  of  E  as  defined 


in  (2.3).  There  are  two  sets  of  unknown  {o..}  and  (6.H  in  E.  A  consistent 

ij  i 


estimate  of  is  provided  by  as  shown  in  Theorem  2.2. 


Theorem  2.2.  Let  F..  be  continuous  at  (0.,0.)  =  (£  .  (q  .  ) , £  .  (q  .  ) ) .  Then 
-  i  J  X  J  MJ 


5..  =  F?n)(c5n)(q.),£(.n)(q.))  =  F^n)(0.,0.) 

ij  ij  Mi  ij  i*  J 


o ,  .  =  F..(0.,0.)  a.e.  as  n  . 
ij  iJ  i  J 


Proof 


| f . . (0 . , 0 . )  -  F!n)(e.,0.)| 
1  IJ  1  J  IJ  1  J 


!  i  Fij(ei»e  j )  -  vvy  I 


-.(n) . 


+  sup|F  (x,y)  -  Fj . ;(x,y) 
x,y  J  J 


Since  F,.  is  continuous  at  (0.,0.)  and 
ij  i  J 


sup|Fij(x,y)  -  F^? }(x,y) 


a.e. 


(2.7) 


(2.8) 


(2.9) 


it  follows  that  the  expression  on  the  left  hand  side  of  (2.8)  -*■  0  a.e.  which 
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establishes  the  result  (2.7)  of  Theorem  2.2. 


The  result  (2.7)  implies  that  o„  in  (2.3)  can  be  consistently  estimated 


by  its  sample  equivalent  . 


There  exist  several  methods  for  the  estimation  of  (see  Kriogor  and 


Pickands  III  (1981)  and  the  references  there  in).  Recently,  a  consistent  and 

-1 


efficient  estimator  of  6^  based  on  a  sample  of  size  n  has  been  proposed  by 
Babu  (1986a)  under  the  assumption  that  f .  is  continuously  differentiable  at 
£^(q^).  There  is  a  possibility  of  this  estimate  taking  negative  values,  and 
when  this  happens  some  modification  of  the  estimate  may  have  to  be  made.  Usii 


consistent  estimates  of  a.,  and  6.  ,a  consistent  estimate  of  o.. /<$.<$.,  the 

ij  i  ij  i  j 


(i,j)-th  element  of  Z,  can  be  obtained  as  0../6.6. 

ij  i  J 


Another  possibility  is  to  obtain  a  direct  estimate  of  o  ../<$.  S  .  by  the 

V1i.) 


bootstrap  method 


S.  ./6.6.  «  E  (n  (9  .  -  0  .  )  (  0  ,  -  0  . )  J 
ij  i  J  1  1  3  J 


(2.10) 


where  E  is  the  expectation  under  the  bootstrap  distribution  function.  The 
consistency  of  the  estimator  (2.10)  can  be  proved  on  the  same  lines  as  those 
given  by  Babu  (1986b)  for  the  bootstrap  estimate  of  the  variance  of  the  sample 
median. 


3.  TESTS  OF  SIGNIFICANCE  BASED  ON  MEDIANS 


Let 


9i 


=  (0 


1  I 


,0  .  ), 

I" 


(3.1  ) 


be  the  marginal  sample  medians  and  an  estimate  of  Z  (as  defined  in  (2.3)) 
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obtained  from  a  sample  of  size  n^  from  a  p-variate  population  IK,  i=l,...,k. 
Further  let  0^  =  (0^,...,0  ^)  be  the  true  value  of  the  marginal  medians  for 
IK.  To  test  the  hypothesis 


8  =...=§ 

1  P 


(3.2) 


we  can  use  the  statistic 


2  k  _ i _  _ ,  ^  _i__ 

X  =  trace[  Z  n.I.  0.0. -(  Z  n.Z.  )09 ' 1 
i=1  i  i  -i-i  i=1  ii  __ 


(3.3) 


where 


k  _  i  ~ 

0  =  (  Z  n.Z.  )  Z  n.Z.  e. 
i=l  i=l  1  1  -1 


(3.4) 


as  chi-square  on  p(k-l)  degrees  of  freedom,  provided  the  individual  sample 
sizes  nj,...,n^  are  large. 

In  cases  where  a  common  Z  for  the  k  populations  can  be  assumed,  we  have 
the  problem  of  estimating  Z  from  the  combined  sample.  For  this  purpose  we 
consider  the  residual  vectors  by  replacing  each  observed  vector  by  its 
difference  from  the  sample  median  vector  computed  from  the  sample  to  which 
the  observed  vector  belongs.  There  are  altogether  n  =  (n^+...+n^)  residual 
vectors,  arising  out  of  the  k  different  samples,  from  which  we  construct  a 
p-dimensional  empirical  distribution  function  E  with  the  marginal  medians  as 
zeros.  Then  o_  can  be  estimated  from  E„,  the  (i,j)-th  marginal  d.f.  of  E  as 
indicated  in  (2.7)  and  <K  from  E^,  the  i-th  marginal  d.f.  of  E  using  any  of 
the  methods  described  at  the  end  of  Section  2.  If  we  denote  a  common  estimate 
of  Z  by  Z,  then  we  can  develop  tests  of  significance  concerning  the  structure 


of  the  median  vectors  0  ,  i=l , . . . ,k,  as  in  the  case  of  mean  values  (see 
Rao  (1973),  p.  558).  For  this  purpose  we  compute  the  "between  populations" 


I 


m 

» 


I 


matrix 


S  =  T.  n  .  9 . 0  .  -  n8G 
111  _ 

i=l  -1-1 


(1.5) 


where  n8  =  n^8j  + 


.  4-  n,8,  ,  and  set  up  the  determinenta 1  equation 

K~  K 


s  -  \i\  =  o. 


(3.6) 


The  roots  of  the  equation  (3.7)  can  be  used  as  in  the  table  on  p.  558  of 


Rao  (1973)  to  test  the  dimensionality  of  the  configuration  of  median  values. 


A.  JOINT  DISTRIBUTION  OF  THE  MARGINAL  QUANTILE  PROCESSES 


In  Section  2  of  the  paper,  we  derived  the  joint  asymptotic  distribution 


of  specified  marginal  quantiles.  We  now  derive  the  weak  limits  of  the  entire 


marginal  quantile  processes  after  suitable  scaling.  More  specifically  we 


consider  the  processes  {Z^}  indexed  by  (q^ . .  )6(0,1)*5  where 


Zn(q1,...,q  )  =  /n[f1(C1(q1))(5jJ(q1)-51(q1)), • • •  , 


f  (ZMjnC' '^(q  )-£  (q  ))]• 
P  P  P  P  P  P  P 


( A  .  1 ) 


We  first  simplify  the  problem  using  the  following  result  which  is  essentially 


a  restatement  of  Theorem  5.2.2  of  CsorgB  and  Revesz  (1981). 


Theorem  3.1.  Suppose  that  for  i=l,...,p,  the  marginal  d.f.  F  is  twice 


differentiable  on  (a^,b^),  where 


*  V.VjXVV 
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-  M  <  a.  =  supfx : F . ( x)=0} 
—  l  1 


00  >_  =  inf  (x  :  ( x  )  =  1  } 


and  F ,  = 


f^  0  on  (a^b^).  Further  assume  that 


max  sup  F. (x)[ 1-F. (x) ]— - 

i  a  <x<b.  1  1  f.  (x) 

i  l  l 


is  non-decreasing  (non-increasing)  on  an  interval  to  the  right  of 


a^(to  the  left  of  b^).  Let 


C(qi . v =  ^v(in)^1)-t*i . vpn)(y-y 


,(n)  , 


where  is  the  empirical  d.f.  of  the  uniform  variables 


.  u..  =  F .  (  x  .  .  )  ,  j  =  1 , .  .  . ,  n . 
n  l  1 1  J 


SUp  J fYn(q)  "  Zn(q)l I  ^  0  a>e' 

qe(0,l)P  "  ~  ' 


(4.2) 


Hence  {Y  }  and  (Z  }  have  the  same  limit, 
n  n 


Note  that  the  marginals  of  {Yn}  converge  weakly  to  a  Brownian  bridge  on 


C[0,1]  (see  Brillingsley  (19G8)).  Since  the  paths  of  the  limiting  process 


are  continuous,  we  define  a  new  process  Y  close  to  Y  as  follows.  Let  D,  (t) 

n  n  1 


be,  as  a  function  of  te(0,l],  the  d.f.  corresponding  to  a  uniform  distributic 


of  mass  (n+1)  over  each  of  the  (n+1)  intervals  (d^  ^ , d ^  J ,  j  =  l,...,n  +  1 


where  d„  =  0,  d  .,=1  and  d,,...,d  are  the  values  of  u.  ,,..., u.  arranged  in 
U  n+1  in  ll  m 


'  _  *■  *  »  .  '  .  •  k'<  . ’i* 


y.v 
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increasing  order.  Clearly 


So  if 


|vJn)(t)-D5n)|  <  -,0  <  t  <  1  a.e. 
i  i  i  '  —  n  —  _ 


Yn(q)  =  ^(Djn)(qi)-qi . Dpn)(^p)'clp) 


P  P  P 


||Yn(q)-Y”(q)| |  <  n"*  Vq e[ 0 , 1 ] P  a.e. 


As  a  consequence, {Y^}  and  {Z^}  have  the  same  weak  limits  and  the  marginals  of 

Y  are  continuous  functions.  Note  that 
n 

Ynt;-B  =  { h :  h  ( q )  =  (h1(qJ  ) . h  (q 

is  a  continuous  function  on  [0, 1  ] ,  i=*=l , .  .  .  ,p} . 

Clearly  B  is  a  separable  closed  linear  subspace  of  the  Banach  space  C^  of 
continuous  functions  on  ( 0 , 1 j P  into  !F^. 

We  shall  show  that  {Y^}  converges  weakly  to  a  Gaussian  measure  on  B.  A 
probability  measure  p  on  B  is  called  Gaussian  if  for  every  H6B  ,  the  space 
of  real  continuous  linear  functionals  on  B,  pH  *  is  Gaussian  on  the  line 
(see  Aranjo  and  Gine  (1980),  pp.  140-142,  28  and  problem  2  on  p.  33). 

'ft 

To  characterize  B  ,  let  H  be  a  real  continuous  linear  functional  on  B. 


H(  h.  ,  .  .  .  ,  h  )  =  H(h.  ,0 . 0)  +  ...  +  11(0,0 . h  ) 

1  p  1  p 


=  H.(h.)  +  ...  +  H  (h  ),  say. 
11  P  P 


(4.3) 
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The  zeroes  in  tie  first  line  of  (4.3)  refer  to  the  zero  function.  Clearly, 
each  is  a  real  continuous  linear  functional  on  C [ 0 , 1  ] .  It  then  follows 
that  B  is  the  k-fold  direct  sum  of  the  dual  space  C  of  C[0,l].  By  Riesz's 
representation  theorem,  for  any  L  e  C  ,  tliere  exists  a  signed  measure  v  on 


[0,1]  such  tht 


:f)  =  [  f ( x  )dv( x ) 
J  n 


for  any  f  e  C [ 0 , 1 ] ,  (see  Dunford  and  Schwartz  (1958)).  Thus  for  every  H  6  B  , 
there  exist  signed  measures  v,,...,v  on  [0,1]  such  that  for  f  =  (f,,...,f  )eB, 

1  p  1  p 


H(f)  = 


P  /-l 

T  f.(x 

i=lJ„ 


)dvi(x) . 


Now  let 


A  =  {  l  a.E  :0<x.£l,x.,a.  rational ,  j=l , . . . , r, . 
j=l  J  xj  J  J  J 

r  =  1,2,...  } 


where  e^  is  the  probability  measure  putting  all  its  mass  at  x.  It  is  easily 
seen  that  A  is  dense  in  C  and  is  countable.  We  now  state  the  main  result. 

Theorem  3.2.  {Y^}  converges  weakly  to  a  Gaussian  random  element 

W  =  (Wp...,W^)  in  B,  where  is  a  Brownian  Bridge  for  each  i  and 


E(W^(t)Wj(s) )  =  P(Fi(xij)<t,  Fj(Xjj)<s)  -  ts 


(4.4) 


for  all  i,j  and  0  <  t,  s  <_  1. 

Proof .  Since  ('/n(D^n  ^(  t  )-t )  :0<t<  1 }  is  tight  for  each  i  in  C(0,1],  it 
follows  that  {Y^}  is  tight  in  B.  Since  A  is  dense  in  C  ,  in  order  to  show  that 


{Y^}  has  a  weak  limit  it  is  enough  to  show  that  for  any  , 

q  q  in  [0,1]  and  a.,  real 

pi  Mpr  ij 


P  r 

£  £ 


i=l  j=l 


a.  ,/n(D^n\v.  .)-q.  .) 
i  ij  ij 


Converges  weakly.  This  holds  because  of  the  central  limit  theorem  and  the 
fact  that 


sup  |V^n)(t)  -  D^n)(t) |  <  -  a.e. 
o<t<l  1  1 


To  complete  the  proof  it  is  enough  to  show  the  existence  of  W  satisfying 
(4.4). 

Since  {Y^}  is  tight,  there  exists  a  random  element  Y  on  B  and  a  subsequence 
{Yn,}  such  that  Y  ,  converges  weakly  to  Y  *  ( Y^ 1  \  .  . . , Y^ ) .  Further,  from  the 
above  arguments 


P 

£ 

i=l 


r 


£ 

J-l 


a. ,Y(l)(q.  .) 
ij  ij 


and 


P 

£ 

i=l 


r 

£ 

j-l 


have  the  same  distribution  as  that  of  normal  random  variables.  So  it  follows 
that  Y  satisfies  the  properties  of  W  mentioned  in  (4.4)  and  Y  is  Gaussian. 

Thus  Yn  converges  weakly  to  W,  and  in  view  of  Theorem  3.1,  {Z^}  converges  to  W. 
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